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ABSTRACT 

We consider the evaporation of close in planets by the star's intrinsic EUV and X-ray 
radiation. We calculate evaporation rates by solving the hydrodynamical problem for 
planetary evaporation including heating from both X-ray and EUV radiation. We show 
that most close-in planets (a < 0.1 AU) are evaporating hydrodynamically, with the 
evaporation occurring in two distinct regimes: X-ray driven, in which the X-ray heated 
flow contains a sonic point, and EUV driven, in which the X-ray region is entirely 
sub-sonic. The mass- loss rates scale as Lx/a 2 for X-ray driven evaporation, and as 

<&y 2 /a for EUV driven evaporation at early times, with mass-loss rates of order 10 10 — 
10 14 g s _1 . No exact scaling exists for the mass-loss rate with planet mass and planet 
radius, however, in general evaporation proceeds more rapidly for planets with lower 
densities and higher masses. Furthermore, we find that in general the transition from 
X-ray driven to EUV driven evaporation occurs at lower X-ray luminosities for planets 
closer to their parent stars and for planets with lower densities. 

Coupling our evaporation models to the evolution of the high energy radiation - 
which falls with time - we are able to follow the evolution of evaporating planets. We 
find that most planets start off evaporating in the X-ray driven regime, but switch to 
EUV driven once the X-ray luminosity falls below a critical value. The evolution models 
suggest that while 'hot Jupiters' are evaporating, they are not evaporating at a rate 
sufficient to remove the entire gaseous envelope on Gyr time-scales. However, we do 
find that close in Neptune mass planets are more susceptible to complete evaporation 
of their envelopes. Thus we conclude that planetary evaporation is more important for 
lower mass planets, particularly those in the 'hot Ncptune'/'super Earth' regime. 

Key words: planets and satellites: atmospheres planets and satellites: physical 
evolution X-rays: stars ultraviolet: planetary systems ultraviolet: stars 
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1 INTRODUCTION 

The majority of planets detected to date are found to be close 
(< 0.5AU) to their parent star. At these small separations, 
and particularly at early times, they will be strongly 
irradiated. The bolometric flux of the host star may inflate 
the planet's radius, above that which would be expected a t 
larger separations, due to heating (e.g iBaraffe et al1l2008l h 
However it is the high energy radiation (UV and X-rays) 
that will be important for the evolution of the planet's upper 
atmosphere, where the gas temperatures may become close 
to the escape temperature. It is important to understand the 
process of planetary evaporation as it may have important 
evolutionary consequences for close-in planets, and it has 
been suggested that gas giants might entirely loose their 
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gaseous envelo pes through such a p rocess, leaving behind a 
rocky core (e.g. I Jackson et aUl2010h . Observational evidence 
for planetary evaporation exists through the detection of 
extended atmospheres in several stellar lines that give 
planetary radii considerably larger than in the optical and 
infrared. In particular the atmosphere extends beyond the 
Roche-radius, indicating the gas is no longer gravitationally 
bound to the planet. The two best studied examples thus 
far are HD 209458b dVidal-Madiar et all 120031 . 120041 ') and 
HD 188733b jLecavelier des Etangs et alj|2010h ' which both 
have estimated mass-loss rates in the range ~ 10 10-13 g s _1 . 

Currently, theoretical models of planetary evaporation 
are limited and are derived under simplistic assumptions, 
additionally there is little consensus whether it is driven 
by Extreme-UV (EUV) or X-ray heating. There is even 
debate as to whether evaporation is in the hydrodynamic 
limit, or proceeds through ballistic loss of particles (e.g. 
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Jean's escape) where the upper atmosphere is no longer 
pressure dominated. Perhaps the simplest approach is to 
assume that every photon received at the pl anet's surface 
is tu r ned into mass-loss at some efficiency (|Watson et al.l 
Il98ll ; iLecavelier des Etangsl 120071 ') . Such an approach can 
only provide an order of magnitude estimate of the mass-loss 
rate, though it is likely to be fairly accurate in the case of 
'energy-limited' evaporation, by which we mean that the 
do minant energy l oss p rocess is PdV work as discussed 
by IWatson et all (|l98ll ). This may indeed represent an 
observable region of parameter space at late times in 
the c ase of EUV driven evaporation |Murrav-Clav et al.l 
120091 ). Nonetheless such an approach masks many of the 
complexities associated with thermally driven evaporative 
winds. Thermally driven winds are free to absorb energy up 
to their sonic surface, whic h may in prin ciple, be located 
far from the launch point (jParkerl 1 1960h . This results in 
the interception of a much higher fraction of the stellar 
high-energy flux than the simple planetary disk, an effect 
lLammer etafl l|2003l ) discussed by introducing a further 
efficiency factor into the 'energy-limited' formalism in the 
form of an expans ion radius based on the work of 
IWatson etal]|l98ll . The inclusion of the expansion radius 
has been dropped in recent years as detailed models of 
EUV heating in hot Jupiters suggested very little energy 
is absorbed high in the flow, with most energy deposited 
near or at th e base of the flow (e.g. Yelle et al. 2004; Garcia 
Munoz 2007: iLeitzinger et alll201lf ). Although it is unclear 
whether EUV heated flows may absorb significant amount of 
energy high in the flow for planets with different properties 
compared to hot Jupiters, and whether the X-rays can 
provide significant heati ng to the gas at large atmos pheric 
heights, as suggested bv lCecchi-Pestellini et ail l|2006l ). 

fn addition, one must account for the reduction in 
the gravitational binding of the planetary atmosphere 
induced by the ste llar gravitational field, as discussed by 
lErkaev et ail (|2007h . In the energy-limited formalism this is 
accomplished by the introduction of another efficiency factor 
based on the assumption one only need overcome the Roche 
potential rather than the gravitational field of the planet in 
isolatiorQ]. Thus th e energy limited formalism presented by 
lErkaev et aD (|2007h gives rise to a mass-loss rate of the form 

• _ ^HE-Rp - . 

m ~ V 4GM p a 2 K(R Roche /R p ) ' 1 ' 

where rj is the 'efficiency' of the flow, Lhb is the high 
energy luminosity of the central star, Mp, R p & a are the 
planet mass, radius and orbital distance and K(Rn oc h c /Rp) 
accounts for the reduction in the planetary binding energy 
due to the Roche lobe. Eq. [1] has been used for several 
parameter studies of both EUV evapor ation (e.g. Lammer 
et al . 2009) and X-ray eva poration (e.g. iDavis fc Wheatlevl 
120091 ; Ijackson et al.l 12012? ); however, the 'energy limited' 
approach possesses tunable parameters, for which there is 
little reason to expect to be independent of the physical 
properties of the system. As such Eq [T] is of limited use in 
understanding the true details of planetary evaporation over 
a range of parameter space. 



Further, there is no reason that evaporation must 
occur in the energy-limited regime and it may instead be 
that radiative losses dominate the energy budget rather 
than PdV work. In the case of EUV e vaporation at 
high luminosities iMurrav-Clav et al l l|2009l ) showed that 
this results in a qualitatively different scaling with both 
the high-energy luminosity and orbital distance. It is also 
necessary to determine whether the planetary evaporation 
is indeed occurring in the hydrodynamic limit, for which 
one needs full solutions to verify that the flow is pressure 
dominated all the way to the sonic surface. 



Several hydrodynamic studies of planetary evaporation 
have been u nderta ken that go beyond the 'e nergy- limited' 
formalism. lYelk 



IMurrav-Clav et al 



1 Note: later we discuss how the effect of the star's gravity appears 
through the derivative of the potential rather than its absolute 
value 



(120041 ). iTian et all (|2005h and 
20091 ) calculate the h ydrod y namic 
escape for a pure EUV heated model, while lYellel l|2004l) 
additionally solves for the hydrogen and helium str ucture 
of the flow. Of these only IMurrav-Clav et all (|2009l ') solve 
for the EUV radiat i on tra nsfer and heating rat e directly. 
However, both lYelld (|2004h and lTian et"all (|2005l ) use fixed 
he ating efficie ncies (for photo-electric heating in the case 
of lYel lc 2004 a nd t he total heating e fficiency in the case 
of ITian et all 120051 '). iPenz et all (^OOct ) include both X-ray 
and EUV, but they set the X-ray heating rate to a fixed 
efficiency, which is then varied betw een 10 and 60 per cent 
to ob tain various flow solutions. ICecchi-Pestellini et al.l 
(|2006h consider pure X-ray heating of a static plane-parallel 
atmosphere and determine the X-ray driven heating with 
consideration of the hydrogen and helium photochemistry. 
They conclude that at high X-ray luminosities similar to 
those found around young stars (e.g. Lx ~ 10 30 erg s _1 ), 
it may be X-rays that dominate the heating. Conversely 
several hydrodynamical studies of EUV evaporation have 
been successfully applied to the observations of HD 
209458b, correctly reproducing the observed mass-loss 
rates of ~ 10 10 — 10 11 g s , suggesting X-rays may be 
unimportant in this particular case (e.g. Garcia Munoj 
|2007| ; IMurrav-Clav et al.ll2009l : iKoskinen etal]|201Cl ) . What 
is not clear, however, is the structure of X-ray driven flows, 
the interactions that take place between a hydrodynamic 
flow that has X-ray and EUV heated regions and how a flow 
may transitions from X-ray dominated to EUV dominated. 



As a prelude to performing multi-dimensional 
hydrodynamic calculations, in this work we build ID, 
on-axis, hydrodynamical solutions to the problem of 
planetary evaporation including both EUV and X-ray 
radiation fields under the approximation of radiative 
equilibrium. This allows us to understand whether 
planetary evaporation will be 'energy-limited' or not, 
whether it will occur in the hydrodynamic limit, and what 
the interactions between the EUV and X-ray fields yield 
in terms of evaporation. Our paper is organised as follows: 
in Section [5] we discuss the various heating mechanisms; in 
Section [3] we derive our hydrodynamic solutions for EUV 
and X-ray heating; in Section|4]we consider simple evolution 
of planets undergoing evaporation; in Section [5] we discuss 
the implications of our results and our assumptions, and 
finally we present our conclusions in Section [6] 
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2 HEATING MECHANISM 

Solar-type stars are known to produce large amounts 
of high-energy radiation in the form of UV and X-ray 
photons. These high-energy photons originate in the stellar 
corona, and at early times may have a luminosity ~ 
1CT 3 of the bolometric luminosity (e.g. iGiidell |2004 ). 
Furtherm ore, this high-energ y radiation falls off with 
age (e.g. IJackson et al.l [20121) in as s ociation with stellar 



spin-down (e.g. Hempelmann et all 1 19951 ; IGiidell |2004| ; 
I Wright et al]|2011 ). Therefore, the impact the high-energy 



photons have on the planetary atmosphere will also vary 
with time. 

Young pre-main sequenc e stars also emit large 
Far-UV (FUV) fluxes (e.g. Ilnglebv et all l201lh . that 
could be re sponsible for drivi ng efficient evaporation from 
disks (e.g. iGorti et all 1201 ll ). However, the main FUV 
heating mechanism to reach the high temperatures (> 
1000K) necessary for planetary evaporation is photoelectric 
heating from dust grains (e.g. iTielens fc Hollenbachl [l~985l ; 
IGorti fc Hollenbachll2003) . Due to the high dust-to-gas mass 
ratios required for efficient FUV heating we will at this 
stage neglect the influence of FUV radiation, though we 
note that FUV heating at lower temperatures may have 
some influence through other heating channels such as H2 
pumping. Thus while it is unlikely that the FUV can be a 
dominant driving mechanism it may induce slight differences 
in the pure X-ray/EUV evaporation discussed here. 



2.1 X-rays 

X-ray photons span a large range in energy hv ~ 0.1 — 
10 keV, with the 'soft' X-rays (hv < 1 — 2 keV) responsible 
for heating in the flux range expected for close-in planets. 
ICecchi-Pestellini et alj (|2006l ) considered X-ray heating for 
a plane-parallel planetary atmosphere with similar fluxes 
to those found around close-in planets, and concluded that 
photo-electric heating could heat planetary atmospheres to 
sufficiently high temperatures to allow evaporation. For 
observed X-ray luminosities, and gas densities in planetary 
atmospheres, heating typically occurs to neutral column 
densities of N = 10 22 cm~ 2 , with the gas reaching 
temp eratures in the r a nge o f a few hundred K to ~ 10 4 K 
Glassgold et all 12004 ICecchi-Pestellini et all l2006l ; 



(c 



lErcolano et al. l2009h . The dominant heating mechanism is 



photo-electron generation by the K-shells of metals, with 
Oxygen and Carbon being the most important. Furthermore, 
the heating time-scale will always be the quickest time-scale 
in the problem (simply the thermalisation time-scale of the 
hot, ~ 10 6 K, photoelectrons) . 

Cooling in X-ray heated gas can occur through several 
processes, though for the expected range of temperatures 
here line cooling will dominate, as X-ray recombination 
and thermal bremsstrahlung are only important at much 
higher luminosities. Line cooling is also dominated by metals, 
and again Carbon and Oxygen are the most important. 
In radiative equilibrium the temperature of X-ray heated 
gas can be described as a monotonic function of the 



Lx/r 



(e.glTarter et al.lll969l ; 



ionization parameter - 

llgea fc Glassgold! 1 19991 ; bwen et al.ll2010h , where L x is the 
X-ray luminosity, n the number density and a the distance 
to the X-ray source. Such temperature-ionization parameter 




Ionization Parameter [erg cm s ] 



Figure 1 ■ Tem perature-ionization parameter relation from 
lOwen et all J2010|V 



(T = /(C)) relations were calculated in lOwen et al.l |2010l) 
for use in discussion of disk photoevaporation around 
solar-type stars. Given the similar ionization parameters 
expected, since the escape temperature from a planetary 
atmosphere is similar to the escape temperature from the 
inner disk (1-10 AU) around a solar-type star, we will 
m ake use of the tem perature-ionization parameter relation 
of lOwen et all l|2010h in this work, as shown in Fig. [T] 

Furthermore, given that the same metals dominate both 
heating and cooling of the gas we may include a basic 
metallicity scaling. To do this we make use of the results of 
lErcolano fc Clarke] (|2010T ) for the heating of protoplanetary 
disks. Interpreting their calculated scaling of the mass-loss 
rate with metallicity (M w oc Z~ ' 77 ) within the theoretical 
framework presented in lOwen et al.l (|2012t ). it is easy to 
show that the ionization parameter scales with metallicity 
as £ oc Z~ ' 77 . Thus, for the simple evaporation models we 
are developing in this work, we incorporate metallicity effects 
using the above scaling. 



2.2 Extreme-UV 

Extreme-UV (EUV) photons have sufficient energy to 
directly ionize hydrogen atoms (hv > 13.6 eV), resulting 
in fully ionized gas at ~ 10 4 K. In radiative equilibrium the 
the ionized gas readily produces further EUV photons by 
recombination, so the gas thermostats to 10 4 K, even in the 
presence of mild heating/cooling from other radiation (e.g. 
X-rays) , thus we use an isothermal equation of state for EUV 
heated gas. The extent of the EUV heated region can then 
be f ound by balancing recomb ination with ionization. 

iMurrav-Clav et al.1 (|2009T ) determined that at low EUV 
luminosities typical of stars with Gyr ages, radiative 
cooling is no longer the dominant cooling mechanism in 
an evaporative flow, and instead cooling through PdV 
work dominates (as much as 80% of the energy went into 
PdV work). At this stage we are interested at looking at 
times when EUV evaporation can contribute significantly to 
the mass- loss from the planet, which will happen at early 
times when the EUV luminosity is high enough to be in 
radiative equilibrium. However, we do take into account 
this 'energy-limited' EUV evaporation when we consider the 
evolution of example planets in Section [4] 
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3 BASIC HYDRODYNAMICS 

In this section we discuss the basics of an evaporative 
hydrodynamic wind from a planet. To do this we 
solve the steady-state hydrodynamic problem along a 
streamline connecting planet and star in a frame co-rotating 
with the planet using simil ar model setups t o previous 
hydrodynamic models (e.g. Tian et al.l 120051 ; iPenz et al.l 



120081 ; iMurrav-Clav et al]|200g ). although we do not assume 
a fixed heating efficiency for the X-rays a re rather calculated 
the gas temperature directly. As in IMurrav-Clav et"aH 



(|2009l ) we neglect the Coriolis force, which can be considered 
as small in the trans-sonic region of the flow as |fl||u| <C 
gravity /pressure. Therefore, the governing equations are the 
steady-state continuity and momentum equations 

'pur 2 ) = 0, (2) 

^du 1 dP G>f cff 

p dr 



Or 



dr 



dr 



(3) 



where r is the distance from the centre of the planet, p is the 
density of the flow, u is the flow velocity, P is the pressure 
and ^cff is the effective (Roche) potential in the co-rotating 
frame, given by 

GM„ GM, 



*„ff = — 



1 G(M„ + M.) 
2 a 3 



M» 



M» + M„ 



(4) 



We adopt radiative equilibrium, which holds in the 
limit of large fluxes (most of the received energy is emitted 
as cooling radiation rather than lost to PdV work), an 
assumption we will justify in detail in Section R^l In a region 
that is exposed to both EUV and X-ray radiation, we adopt 
an isothermal equation of state as discussed in Section 12.21 
given by 

Peuv = pc| uv , (5) 

where ceuv is the isothermal sound speed in ~ 10 4 K gas, 
for which we adopt a value of 10 km s _1 . In gas that is 
only exposed to X-rays, we adopt a temperature-ionization 
parameter relation, T = /(C) 1 where / is a monotonic 
function taken from t he previous X-ray calculations of 
lOwen et al] (|2010l . l201lh , as disused in Section [2~lffl . 

The transition from an EUV to X-ray heated wind 
will occur at an ionization front located at a distance 
r = -Rif from the planet. S uch a flow is q u alitat ively 
similar to those described by Ijohnstone et al] (| 1998T ) for 
the photoevaporatio n of disks illuminated b y external UV 
radiation, where in Ijohnstone et al.l (|l998l ) an ionization 
front separates an EUV heated flow from an FUV heated 
flow. The location of the ionization front relative to 
the sonic point is crucial in determinin g which heatin; 
mechanism sets the mass-lo s s rate (e.g 



Johnstone et al 



f 



Il998l : iHollenbach et al 1 l2000l : lOwen "et al.ll2010h . since the 
flow above the sonic point is causally disconnected from that 
below it. Therefore, in a similar manner to the FUV or E UV 
dominated photoevaporation of Ijohnstone et al] ljl998h . wo 
can identify two separate cases of planetary evaporation: 



2 We note that this function is empirical and simply use the label 
/ to denote it in this work 



X-ray dominated, where the sonic transition occurs in the 
X-ray heated flow, and EUV dominated, where the ionization 
front occurs before any sonic transition in the X-ray heated 
region of the flow. In the latter case a partially neutral, 
sub-sonic, X-ray flow, passes through an ionization front 
into an ionized, transonic, EUV flow. These two cases are 
illustrated schematically in Fig. [21 where the left hand panel 
shows X-ray driven evaporation and the right hand panel 
EUV driven evaporation, where the EUV heated region is 
always closer to the central star as the X-rays have a much 
lower absorption cross-section. 

Here, we are only interested in mass-loss rates, which 
do not depend on the structure of the flow beyond the sonic 
point, so we do not solve for the EUV/X-ray transition in the 
supersonic region of an X-ray dominated flow. In reality the 
X-ray wind may shock and form another sub-sonic X-ray 
flow, which then passes through an ionization front into a 
supersonic EUV heated flow, or it may just supersonically 
pass through a rarefaction ionization front directly into 
a supersonic EUV heated flow. The exact configuration 
(shock or no shock) will be deter mined by the M ach number 
upstream of the ionization front (|Spitzerlll978l ). In an EUV 
dominated flow, however, we do solve for the EUV/X-ray 
transition, as the X-ray heated flow will be entirely sub-sonic, 
and the associated possible sonic point will be in the EUV 
heated region, determined by the escape temperature at the 
ionization front. The position of the ionization front is set 
by balancing recombination and ionization at the ionization 
front with the mass-flux (rhx) passing through the ionization 
front. That is 

2 



<I>, 



Ana 2 



T6tt 2 



uceuv 



(6) 



where <E>* is the luminosity of EUV photons (photons per 
second), a r = 2.6 x 10~ 13 cm 3 s _1 is the recombination 
coefficient for hydrogen at 10 4 K, p is the mean particle 
weight, and A is an order unity geometry factor that 
takes into account how steeply the density falls off in the 
ionized portion of th e flow, with typically A ~ 1/3 (e.g 
Ijohnstone et al] 1 19981 ). Therefore, since Eq. [6] gives _Rif oc 
m^/ 3 , there is a critical rhx a t which the ionization front 
occurs at the same radius as the X-ray sonic surface (Rw = 
R a ). At an rhx below this critical value the flow will be 
EUV dominated, as the X-ray portion of the flow is entirely 
sub-sonic, whereas at an rhx above this critical value the flow 
will be X-ray dominated, as the sonic point occurs before 
Rif. As we shall show in Section l3.il the sonic surface for an 
X-ray dominated flow typically occurs within a few planetary 
radii of the planet, so if we characterise R 3 by R a = /3R P , 
then /3 is of order unity. Putting astrophysical numbers in 
Eq. [6] we can determine the EUV luminosity necessary for 
the ionization front to occur before any sonic transition in 
the X-ray flow, and thus for which the flow will become EUV 
dominated as 



> 10 40 s" 1 



0.1 AU 



mx 



(— 

10 12 g s- 1 J \l/3 



Rp 



1.5 J \10Rfi 



(7) 



Given these values we can expect that there will be regions 
of parameter space corresponding to either EUV driven flows 
or X-ray driven flows as both can occur under sensible 
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Possible Shock 



EUV 



X-ray / X-ray / 




Ionization Front Sonic Point 

(a) X-ray driven evaporation 





Possible 
Sonic Point Ionization Front 

(b) EUV driven evaporation 



Figure 2. Flow topologies for the X-ray driven and EUV driven evaporation cases. In the case of X-ray driven evaporation the X-ray 
heated flow passes through a sonic surface, then possibly shocks before passing through the ionization front and becoming EUV heated. 
In the case of EUV driven evaporation a sub-sonic X-ray flow passes through the ionization front, the EUV heated flow then is either 
supersonic or proceeds to pass through an EUV heated sonic point. See text for a discussion of whether the shock exists in the case of 
X-ray driven evaporation, and whether the sonic point exists in the case of EUV driven evaporation. 



physical conditions. As we shall show later rhx oc Lx/a 2 

3/2 

and will scale with planet radius steeper than R p , thus we 
can expect EUV dominated evaporation to occur for small 
planets, low X-ray luminosities and at large distances from 
the central star. 



3.1 X-ray driven 

If we assume the sonic point occurs before the ionization 
front, then we can proceed in calculating the transonic flow 
solution. As mentioned above, at this stage we are purely 
interested in obtaining mass-loss rates, so we ignore any 
interactions with the EUV in the super-sonic portions of the 
flow (we are essentially setting $„ = in this sub-section). 
By making use of the fact £ = Lxp/pa. 2 , and given a is 
fixed for planets in approximately circular orbits, we see 
immediately that temperature is purely a function of density, 
or more precisely the flow is barotropic0. 



3.1.1 Basic Structure 

Before using the reali stic temperature-io nization parameter 
function calculated by I Owen et al.ll|2010h . it is instructive to 
first set /(£) oc £ a in order to analyse the basic properties of 
an X-ray driven flow. Using Eq. [2] and the ideal gas equation 
we can re- write Eq. [3] in the form 



(1 - a) c 2 x ] 
= (1 



<91oe 



a) c x 



<9* cff 



(1 



Or 



(8) 



where we define c 2 x as an isothermal sound speed for a gas 
temperature T. Eq. [5] indicates that we can only obtain 
accelerating solutions (du/dr > 0) if a < 1. This is trivial to 
understand, as one is unable to obtain a negative pressure 
gradient in an outflow if a ^ 1 (91ogP/91ogp < 0). For 
a < 1, Eq. [5] possesses a critical point at u = \j 1 — a cx , 



3 The gas behaves as if it has a ratio of specific heats, 
dlog/ 



which can be identified as the sonic point, such that to obtain 
a sub-sonic to super-sonic transition one requires that 



2 

c x = 



1 Rs d'fcB 



dr 



(9) 



Such a condition has a very important effect - it fixes the 
temperature, and hence density, at the sonic surface. In 
essence, since p and u axe specified at this fixed radius, the 
mass-flux is also fixed here. 

Eq. [9] also indicates that there is a maximum radius 
at which the sonic surface may occur, as one requires 
d^eg/dr > 0. The Roche-radius thus ultimately represents 
the maximum radius for a thermally driven outflow. 
Obviously planets with radii greater than the Roche-radius 
will be undergoin g dynamical Roche-lobe overflow (e.g. 
Ijaritz et all 120051 ). rather than wind driven mass-loss. 
Furthermore, we can now calculate the entire flow solution 
by making use of the fact the Bernoulli function is conserved 
along a streamline in a barotropic flow (so numerical 
integrations of the fluid equations are not required). Solving 
for the stream-function of a streamline connecting the planet 
and star we find 



u 

T 



lR^ 9*eff 

q 2 dr 
1 

2 ~~ 



— +*eff(r) 

r=R s \PsJ 

2 dr r=Rs 



(10) 



+ ^ eff (R a ), 



where the RHS of Eq. [10] has been evaluated at the sonic 
surface. Eq. [TO] specifies a family of flow solutions which 
depend on the chosen location of the sonic-surface, R s . In 
reality the chosen flow solution, and hence chosen R 3 , is 
set by dynamical balance in the upper atmosphere of the 
planet, and where the atmosphere transitions from being 
X-ray heated to heated by the star's bolometric luminosity 
(which we identify as the planetary radius, R p ). In Fig.[3]we 
show density, velocity, pressure and column density profiles 
for the solution to this simple problem, varying the planetary 
density, mass, and the index a. From this we can see that 
lower values of a result in broader profiles, as do lower 
planetary masses and densities. 

One of the important features of flows with an effective 
ratio of heat capacities < 1 is the steep gradients in the very 
sub-sonic portions of the flow, close to the planet. In these 
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Figure 3. Density, velocity, pressure and column density profiles for solutions to the X-ray evaporation problem, for various values of 
the temperature ionization-parameter slope (a=0.2, 0.4, 0.6 & 0.8). All profiles are scaled to their value at the sonic surface. The solid 
line indicates the sub-sonic portion of the flow while the dashed line shows the super-sonic portion. The first column shows the solutions 
for a Jupiter mass planet with density 0.5 g cm -3 , the second column a Jupiter mass planet with density 1 g cm -3 , and the final column 
a Neptune mass planet with density 0.5 g cm -3 . 



very sub-sonic regions we can solve for the density structure 
by neglecting the kinetic energy term in the flow (first term 
on the LHS of Eg. I10p . and assuming the effective potential 
is given by $ = —GM p /r. Solving for the density structure 
in the sub-sonic flow in this way we obtain 



, 3 
l + -a 



-l/a 



(11) 



The pressure scale height of this atmosphere can be 
compared to the scale height of a typical underlying 
bolometrically heated planetary atmosphere. Expanding 



about R p it is easy to show 



dlogP 
d log r 



<r-R P 



(12) 



Since this diverges as one approaches R p , it will always tend 
to a value considerably steeper than the typical scale height 
of the planetary atmosphere (> 100km). This implies that 
while the position of the sonic surface will be sensitive to the 
planetary radius, the flow solution will be highly insensitive 
to the structure of the planetary atmosphere below. For a 
given change in the penetration depth of the flow (e.g. due 
to a change in the planet's atmospheric structure, or the 
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X-ray luminosity) APfl ow 3> AP atm , so a large change in the 
structure of the planetary atmosphere causes a very small 
change in resultant flow solution. Thus, the evaporation rates 
are sensitive to the intrinsic properties of planetary mass 
and radius, but very insensitive to the atmospheric structure 
at R p . As such, though the material composition of the 
atmosphere may be important through metallicity effects in 
the flow, we can neglect any details of the structure of the 
planetary atmosphere below the flow. 

With this in mind, and taking R p to be the point at 
which density diverges towards infinity, we can estimate 
Rs/Rp in terms of a. Setting p — > oo in Eq. [TT] we find 



Rs = 



2a 



Rp- 



(13) 



Thus we see that unless a < 1 the sonic surface generally 
occurs within a few planetary radii, and that X-ray heating 
cannot support flows with sonic surfaces at very large 
distances from the planet surface, unless the flow is close to 
isothermal. This is what will ultimately limit the attainable 
mass-loss rates, and it is unlikely that the planetary wind 
can have a sonic surface close to the Roche-radius. 



3.1.2 Realistic Solutions 

The flow solutions discussed above are qualitatively useful 
in understanding the features of X-ray flow, but should 
not be taken as accurate, representative, models for the 
evaporation of planetary atmospheres. In order to calculate 
more representative solutions, we need to use a realistic 
temperature-ionization parameter relation to calculate the 
flow solutions, as discussed above. In Appendix A we show 
that it is possible to obtain a family of transcendental 
solutions for the ionization parameter in the flow for different 
choices of the sonic point, using the same techniques used to 
obtain Eq. [10] 

To determine the position of the sonic surface for a given 
planetary radius we can take advantage of the insensitivity of 
the flow solutions to the details of the underlying (non-X-ray 
heated) atmosphere, and instead choose a sonic radius and 
solve for the planetary radius as the location where the 
pressure in the X-ray flow diverges (as in Eq. I13[) . The 
correct solution, and sonic surface position, is thus that 
which reproduces the desired planetary radius. The accuracy 
of this method has been tested against more detailed 
(and computationally expensive) calculations that directly 
account for dynamical balance in the upper atmosphere 
of the planet. We find that the error in the sonic radius 
position, and mass-loss rate, is small (< 1%) over large, order 
of magnitude changes in both temperature and planetary 
density, covering a range of parameter space larger than 
expected for exoplanets. In addition we find that the sonic 
surface always occurs within optically thin, X-ray heated, gas 
(N(R S ) < 10 22 cm -2 ), validating the ionization parameter 
method. 

Examples of such solutions are shown in FigureUfor two 
planetary radii (« 1.3 & 3.5 x 10 10 cm/ 19.8 & 55 R®) for a 
1 Jupiter mass planet at a separation of 0.1 AU from a solar 
mass star with Lx = 10 30 erg s _1 , although all flow solutions 
show very similar features. We again note the steep gradients 
close to the planet surface, and the sonic point occurs at 




0.5 1 1.5 2 

log 10 (Planetary Radius) [R ffi ] 

Figure 5. Mass-loss rates for planets with solar metallicity, 
located at 0.1 AU around a solar mass star with Lx = 
10 30 erg s _1 . The thick solid line shows the maximum radius at 
which a flow can occur (Roche radius), with the grey shaded region 
corresponding to dynamical Roche-lobe overflow. The magenta 
dashed lines are lines of constant average density; 0.2, 1, 5 g cm -3 
(Jupiter and Earth have densities 1.3 g cm~ 3 and 5.5 g cm— 3 
respectively), and the black dotted lines are 'evaporation lines' 
for ages of 100 Myr and 1 Gyr. The black squares show a sample 
of currently observed transiting planets as a guide to the range of 
masses and radii observed for planets. 



only a few planetary radii in both cases, while at large 
radii the flow is approximately isother mal and share s many 
properties with Parker's thermal wind (|Parkerlll95Sh . Given 
that the flow solutions are specified uniquely as a function 
of the ionization parameter, £, along the streamline, and the 
flows are highly insensitive to the atmospheric structure, 
we can use similarity arguments to extract exact scaling 
relations for several of the flow variables. Since £ oc Lx / p and 
solutions with the same £(r) are topologically identical for 
fixed parameters - M p , R p , a & M* - then £(r) is fixed along 
the streamlines if p oc Lx ■ In the case where the sonic surface 
occurs far from the Roche radius, which is almost always the 
case for realistic planetary densities (so the effect of stellar 
gravity and the centrifugal force is small), flow solutions 
with the same £(r) profile are topologically identical for fixed 
parameters (M P ,R P ). Therefore, the density in the flow then 
scales as p oc a~ 2 in order to yield a fixed £(r) profile along 
the streamline. Furthermore, we can include a metallicity 
scaling using the £ oc Z~ ' 77 result mentioned in Section [2] 
Thus we similarly find p oc Z-°- 77 in a flow with fixed £(r). 
Combining these we find a mass-flux that scales as 



rax oc 



L} 



i z o: 



(14) 



provided R p is not close to the Roche radius, however, exact 
scaling relations do not exist for R p and M p and these must 



© 2012 RAS, MNRAS OOO.mflSl 



8 James E. Owen & Alan P. Jackson 




Radius [cm] 



Figure 4. Profiles of the flow structure for an X-ray heated flow using the realistic temperature-ionization parameter relation, for two 
planetary radii of ~ 19.8 & ~ 55 -Rgj (rs 1.3 & 3.5 X 10 10 cm) for a Jupiter mass planet at 0.1 AU around a solar mass star with 
Lx = 1 X 10 30 erg s — 1 . The solid line shows the sub-sonic part of the flow, while the dashed line shows the super-sonic portion. We note 
all flow structures are similar to those shown here. 



be calculated numerically. It is worth pointing out that our 
definition of 'close' to the Roche-radius is somewhat different 
to what is meant in previous energy-limited calculations (e.g. 
lErkaev et al.l [20071 ). as the hydrodynamic solution actually 
depends on the gradient of the effective potential (Eq. A2) 
rather than its absolute value. So by 'close' we mean troche - 
R s < O.IRp, which, given R s is similar to R p and most 
planets do not have radii that fill their Roche-lobes, means 
Roche-lobe effects to the scalings of these hydrodynamic 
solutions will only matter in the most extreme cases. 

Importantly, as mentioned above, we must also check 
that our flow solutions are in the fluid limit. Physically we 
require that the exospheric radius (where the fluid is no 
longer pressure dominated) is at a radius greater than the 
sonic radius, since anything that happens after the sonic 
point is not in causal contact with the sub-sonic portion of 
the flow and cannot affect it. To do this we compare the 
pressure scale height £ = dr/d log P up to the sonic surface 
with the collisional mean free path of the gas particles given 
by A = 1/na. If i/\ > 1 up to the sonic surface then the 
flow will be in the fluid limit, whereas if t/X < 1 then the 
calculated fluid solution is not appropriate and mass-loss is 
likely to occur from a ballistic type escape mechanism (e.g. 
IVolkov et aI1l201lD . 

In Fig. [5] we show a slice through the Mp and Rp 
plane, showing the mass-loss rates as a colour map. The 
thick solid line indicates the Roche radius, the maximum 
radius at which a flow can occur, with any planet beyond this 
in the grey shaded region undergoing dynamical Roche-lobe 
overflow. One consequence of R s being slightly larger than 
R p is that there is a range of planetary radii close to the 
Roche-radius where no transonic solutions exist (white-space 
between the colour map and thick black line), as any sonic 
surface would have to occur outside the Roche-radius, which 
as discussed above is unphysical. In reality we hypothesise 
that a planet in this region may drive a purely sub-sonic 
wind to the Roche-radius and undergo a thermally driven 
Roche-lobe overflow in this way. It is however important to 
realise that this range of parameter space is extremely small. 
The white region at the left of the figure indicates the region 
where hydrodynamic flow solutions do not exist and instead 
mass-loss may occur via a ballistic type mechanism. This 
figure clearly indicates that the most interesting regime for 
planetary evaporation is the Neptune/super-Earth regime, 



and that in some cases planets of Jupiter mass and 
above may not experience hydrodynamic evaporation at a 
separation of 0.1 AU. 

We also show two 'evaporation lines' in Fig. [S] at times 
of 100 Myr and 1 Gyr for which we set rhxt = M p . A 
planet close to these lines can thus be expected to have 
lost a substantial fraction of its mass in 100 Myr or 1 Gyr 
respectively. As a guide we show a sample of close-in planets 
from tak en the Exoplanets E ncyclopedijfl on 13 September 
2011 (see I Jackson etal1l2012l for a more detailed description 
of this sample) . It is important to note though that we have 
not accounted for the true orbital distances of the planets 
here, they are intended purely as a guide for the range 
of masses observed (see Section [4] where we also take into 
account the separations). 

Including all the effects described above we can calculate 
flow solutions for a range of input parameters and we show 
the variation of the resultant mass-loss rates in Fig. [6] where 
we plot colour maps in the R p — M p plane for various 
values of a, Lx , M*, and Z. Fig. [6] shows that, as expected, 
planets of larger radius for the same mass have larger 
photo-evaporation rates, with the largest values occurring for 
planetary radii just inside the Roche-radius. Furthermore, 
at higher X-ray luminosity and smaller separations the 
hydrodynamic limit of evaporation extends to smaller 
planetary radii. In addition, as in Fig. [5] we see that while 
Jupiter mass planets may undergo evaporation for a small 
range of parameter space, Neptune and super-Earth type 
planets with sensible densities can experience a reasonable 
level of hydrodynamic evaporation during their lifetimes. 

We can compare our results here to the previous 
energy-limited evaporation models summarised by Eq. [JJ 
and find that Eq. [1] performs poorly in reproducing our 
results unless the energy efficiency (77) and expansion radius 
(/3) are allowed to vary by large amounts (in some cases 
orders of magnitude) . Additionally, as we discuss in Section[5] 
we find that X-ray evaporation is not actually occurring in 
the energy-limited regime. Thus, while the mass-loss in the 
X-ray driven regime does have the same scaling with X-ray 
luminosity and separation (Eq. I14|l , we find that Eq [JJ is of 
limited use for X-ray driven evaporation, particularly when 
considering wide ranges in planetary mass and radius. 

4 www.exoplanet.eu, see also lSchneider et al. I (l201lh 
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Figure 6. Mass-loss rates in the Mp-Rp plane, varying the input parameters Lx, a,, M*, and Z. In each panel all parameters except 
the one indicated is held fixed at Lx = 10 30 erg s —1 , a = 0.1 AU, M* = IMq and/or Z = IZ@. As in Fig.[5]the thick solid line indicates 
the Roche radius, while the dashed line indicates a mean planetary density of 1 g cm -3 . 



3.2 EUV driven 

In the case of EUV driven evaporation, the ionization 
front is below the X-ray sonic surface and thus any 
X-ray flow is sub-sonic. As discussed previously EUV 
driven evaporation can exist in two limits where ionization 
is balanced by either gas advection ('energy-limited') or 
recombination ('recombination-limited'). At high fluxes, 
similar to the values during the star's early evolution, EUV 
evaporation can be consid ered to be recombination-limited 
(|Murrav-Clav et all 120091 ). In order to assess the impact 
of evaporation on the planet mass, at this early stage we 
are concerned with the high luminosity phase of the stellar 
evolution and consider EUV driven evaporation to be in 
the recombination-limited regime. In this case the ionization 
front is set by balancing ionization against recombination. 



An ionization front can be of two types; when the 
speed on the neutral side V[ n > 2ceuv the ionization 
front is a rarefied (R-type) front, whereas when v ln < 
cx/2(cx/ceuv) the ionization front is a detached (D-type) 
front (Spitzer 1978). In most cases ceuv S> cx, so the 
ionization front will be D-type and there will be a density 
jump across the ionization front with a maximum input 
speed fin, max = cx/2(cx/ceuv). As such the neutral X-ray 
flow will be highly sub-sonic and the density structure will 
be approximately described by the hydrostatic structure. 

As described in Section fe-ll such a hydrostatic structure 
will be very steep near the planet surface. As such we cannot 
assume that the sub-sonic flow occurs over of order one 
scale height or u se a total colum n argu ment to determine 
its width as in Ijohnstone et al.l 1 19981 . Furthermore, the 
speed on the ionized side of the ionization front will 
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depend on escape temperature at the ionization front (see 
Figure 2(b) I. If the escape temperature is less than 10 4 K 
at the ionization front, then there is no sonic point in 
the flow and the speed on the ioniz ed side of the front is 
fout ~ ceuv (jjohnstone et al.1ll998l ). However, should the 
escape temperature at the ionization front be > 10 4 K then 
the speed on the ionized side of the front will be described 
by the Parker wind solution (v ou t = /parkerCEUv); as a result 
of the necessity that an isotherm al flow has a sonic p oint at 
R- 1 = l/(2c E uv) d^ cft /dr\ r=Rs (|Parkerll 19581 . 1 1960 1. 

In the recombination limit we are considering here the 
position of the ionization front is set by balancing ionization 
and recombination along a line from the star to the planet, 
so 



A-KCL 2 



a r n EUV = Aa r n EVV RiF- 



(15) 



Now writing Rif — (1 + x)R p , where x encapsulates the 
stand-off distance of the ionization front due to the partially 
neutral X-ray flow behind it, the total mass-loss rate may 
be written as mEuv = 47r/mEuv/ P arkerCEUv-RiF- Evaluating 
this at the ionization front we obtain the following equation 
for the total mass-loss rate in EUV driven planetary 
photoevaporation, 



m E UV = 2.4 X 10 11 g S 1 (1+X) 2 /parkcr 



10 40 s- 



O y 1 f Rp \ ' 2 f A_\ 2 / CEUV 

0.1AU/ \WR 9 J \1/3J VlOkms" 1 

(16) 

We note that Eq. 1161 appears throughout the literature for 
EUV driven transonic winds. For the case of /parkcr = 1 
and x = this solution is identical to the well known 
solution for the evaporation of gas c louds of radius Rp 
presented by iBertoldi fc McKed (|l990h . This solution also 
appears for e xternally irradiated pro toplanetary disks where 
/parkcr = 1 (jjohnstone et al.l ll998T ) , and for hot Jupiters 
irradiated wi th a high luminosity p ure EUV spectrum 
with x = (|Murrav-Clav et al.l [2009V As can be clearly 
seen, in the case of EUV driven evaporation there is a 
qualitatively different scaling with the stellar high-energy 
luminosity and the orbital distance than for X-ray driven 
evaporation (rhx oc Lx/a 2 ). This indicates that at large 
separations and low luminosities it will be the EUV rather 
than the X-ray that is dominating the evaporation rates. 
To find x we must solve the jump conditions across the 
ionization front. Using mass and momentum conservation, 
and writing the pressure as Px = (?xPx on the X-ray heated 
side and Peuv = CeuvPeuv on the EUV heated side, the 
deterministic equation for the density behind the ionization 
front becomes 



/parkerCEUv ( ) + /(Lx/nxO 2 ) 



V nx 



— Teuv (l + fp, 



n, ( riEuv 

arker 



V nx 



(17) 



where «euv is given by Eq. \T5\ This equation possesses two 
roots, although only one can correspond to an out-flowing, 
sub-sonic X-ray heated flow0. Once nx is found a Bernoulli 

5 There always exists the possibility of an inflowing solution. 



potential can be constructed for the sub-sonic portion of the 
flow, and x can then be found in an identical fashion to 
that used to find the full X-ray solutions in Section 13.11 In 
Figure[7] we show an example solution of an EUV driven flow 
above a sub-sonic X-ray flow from a Jupiter mass planet. 

In order to consider the role the X-ray radiation plays 
in enhancing the mass-loss rate we show both the total 
mass-loss rates and the value of x for a range of parameter 
space in Fig. [8] This figure shows the mass-loss rates and 
stand-off distances from planets located at 0.1 AU around 
a solar mass star with = 10 40 s _1 , where Lx has been 
lowered from the nominal value used earlier in this work to 
10 29 erg s _1 , to increase the parameter space over which the 
EUV dominates the evaporation. We note that in most cases 
i« 1 although x does become of order unity at higher values 
of Lx and larger planetary radii. In this case the white region 
between the colour-map and the Roche radius is in part due 
to EUV driven solutions not being possible in this region as 
the flow has transitioned from EUV to X-ray dominated. 

To explore the parameter space further Fig. [9] shows 
the variation of the mass-loss rates and region of EUV 
dominance with $*, Lx, and a, where deviations from the 
input parameters used in Fig. [8] are shown in the panels. 
This shows that, as expected due to the different scalings of 
the mass- loss rate from EUV (m E uv oc <&V 2 /a) and X-ray 
(riix oc Lx/a 2 ) dominated flows, EUV photo-evaporation 
dominates at low X-ray luminosities and large separations, 
whereas X-ray photo-evaporation dominates at high X-ray 
luminosities and small separations. As a result X-ray 
dominated winds will in general lead to larger mass-loss rates 
than EUV dominated winds. 



3.3 Direct comparison to HD209458b 

HD209458b has an extended atmosphere which has 
been well studied in the pa s t, pa rticularly using pure 
EUV evaporation. IPenz et al.l (|200&f ) estimated an X-ray 
luminosity of ~ 1.1 x 10 27 erg s _1 and an EUV luminosity 
of ~ 6 x 10 37 s _1 (their Lxuv - L x ) for HD209458b, 
such a combination places it firmly in the EUV driven 
regime as an X-ray luminosity of ~ 10 29 erg _1 would be 
required to switch to the X-ray driven regime. Therefore, 
we indicate that the previous pure EUV approach for 
modelling HD209458b's evaporation is probably accurate 
and our models indicate a mass-loss rate of ~ 2 x 10 10 
g s _1 in good agreement with p revious EUV only models 
and observations of HD209458b (|Vidal-Madiar et alJliooll , 
2004). Additionally we find the X-rays only make a small 
difference to the flow structure at small radius with an 
sub-sonic X-ray heated flow of size < 0.01i? p . Given the 
very small correction the inclusion of X-ray heating makes 
when considering the current evaporation of HD209458b 
in this case our conclusions essentially reduce to a pure 
EUV flow similar to those calculated by iMurrav-Clav et al.l 
(|2009h . As a result of the simplifying assumptions made 
in the treatment of the EUV in this work, previous EUV 
only models that include more detailed chemistry (e.g. 
iKoskinen et alj|20ld :l Ben-Jaffel fc Sona Hosseinil2010h may 
better describe the details of the flow solution in this case. 
At early times however, HD209458b will have undergone 
evaporation in the X-ray dominated regime with mass-loss 
rates ~ 10 13 g s _1 , several times larger than the predictions 
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Figure 7. The structure of an EUV driven flow, with a sub-sonic X-ray heated flow at smaller radius. The X-ray heated flow is shown 
by the solid line, the EUV heated flow by the dashed line and the ionization front by the dotted line. The panels show density, velocity 
and temperature (from left to right), as a function of radius. The planet is 1 Mj at a separation of 0.1 AU, with a radius of 12.9 R^ 
(8.2 X 10 9 cm), an X-ray luminosity of Lx = 10 29 erg s — 1 and an EUV luminosity of $* = 10 4() s — 
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Figure 8. Mass-loss rates (left) and stand-off distance of the ionization front (right) in the Mp — Rp plane. As in Fig. \E\ the planet is 
solar mctallicity at 0.1 AU from a solar mass star and in the left hand panel we include the same lines of constant density and sample of 
known transiting exoplanets. The star has an EUV luminosity, $* = 10 40 s — 1 , while the X-ray luminosity is lowered compared to Fig. [5] 
to Lx = 10 29 erg s — 1 to increase the parameter space over which the EUV dominates the evaporation. 



of previous EUV models, indicating that HD209458b may 
have lost up to 2 — 5% of it's original ma ss, again somewhat 
large r than previous EUV only models (|Murrav-Clav et al.l 
2009 found ~ 0.6%), but certainly not enough to have made 
it unstable to complete destruction. 

In the case where detailed models show that the heating 
is EUV dominated and predominately takes place in a small 
layer near the base of the flow, Eqn. [1] may be used to 
provide reasonable estimates of the mass-loss rate with an 



appropriately chosen efficiency (as in the case of HD209458b 
discussed above r\ ~ 25% e.g. Yelle et al. 2004) . However, we 
caution that the heating efficiency certainly does not remain 
constant for the parameter space considered in this work and 
as such care must be taken when using Eqn.Q]to estimate the 
mass- loss rate. Where a choice of an appropriate 'efficiency' 
requires detailed calculations, such as in this work or others 
(e.g. Yelle et al. 2004; Murray-Clay et al. 2009). 
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Figure 9. Variation of the EUV driven evaporation rate with various input parameters. Deviations fron the standard parameters of 
$* = 10 40 s _1 , Lx = 10 29 erg s -1 , a = 0.1 AU and M* = IMq are shown in the panels. Lines as described in Figure[6] White space 
between the colour-map and the Roche lobe is largely due to transition from EUV domination to X-ray domination. 



4 EVOLUTION 

To assess the effect of evaporation on the planet's evolution 
in this section we build simple evolutionary models to 
indicate whether there are any regions of parameter space in 
which evaporation can significantly influence a planet over 
its lifetime. To do this we follow the evolution of the planet 
mass through time, as well as th e evolution of t he ste llar 
high-energy spectrum. We follow Ijackson et al.l (|2012l 'l to 
include the evolution of the stellar high-energy spectrum, 
using their results for the evolution of the X-ray component 
of the spectrum. After an initial phase (t < r sa t) where Lx 
is saturated at a level L s x 1 , the X-ray luminosity then decays 



L x (t)=LT (t/r sa 
where values of L 



(18) 



T sat and the power la w index a for 
different spectral types are listed in Table 3 of ljackson et al.l 



<|2012h . In the case of pure X-ray evaporation we can use 
the results of Section 13.11 to estimate a planetary survival 
line at ~3 Gyr for different separations under the very crude 
assumption of evolution at constant mass and radius (see 
below for models that follow the evolution of the planet 
mass), to compare to observations. This is shown in Fig. 1101 
where models at four separations, 0.0125, 0.025, 0.5 and 0.1 
AU, are shown. Any planet above the corresponding survival 
line should have completely lost is envelope to evaporation 
at an age of ~3 Gyr, while a planet on or just below the 
survival line would indicate that it has experienced a period 
of strong evaporation. Planets well below the survival line 
would indicate that while the planet might be evaporating, 
it has not changed significantly from its initial mass. 



The 
models 

line. As mentioned above the observations also indicate that 



observations indicate that, as expected if our 
are accurate, no planets lie above their survival 
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the most interesting regime for planetary evaporation is 
the hot-Neptune regime, and that the hot-Jupiters, while 
evaporating, do not loose substantial fractions of their initial 
mass. 



4.1 Transition from X-ray to EUV driven 

The evolution of the EUV luminosity is still observationally 
rather unclear. Motivated by the fact that EUV luminosities 
and X-ray luminosities seem to have similar values at earlier 
times during the T-Tauri phase (lAlexander et all 12004 ) 
and at late times (|Ribas et al.l 120051 ). perhaps - indicates 
fall off with time is sim ilar, if not identical, (see also 
ISanz-Forcada et al.l 1201 ll ). we make a simple assumption 
and couple the X-ray and EUV luminosity to have identical 
values at all times. Thus we take the EUV luminosity to fall 
in step with the X-ray luminosity. 

The transition from X-ray driven to EUV driven 
evaporation is both of observational and theoretical interest. 
Using the models calculated in Section [3] we can estimate 
the X-ray luminosity at which a planet of given mass, radius 
and separation will switch from becoming X-ray driven to 
EUV driven. The results of these calculations are shown in 
Fig. [TT] We note that the X-ray luminosity at which this 
transition takes place for Jupiter mass planets with densities 
of ~1 g cm -3 at separations of > 0.05 AU are considerably 
higher than any expected X-ray luminosity for a solar-type 
star and are not shown. As expected, at higher planetary 
densities, and larger separations, the transition takes place 
at higher X-ray luminosities. In Section f4. 21 we will see that 
this has implications for whether the EUV or X-ray will drive 
the evaporation at late (~Gyr) times, at which we observe 
most planets. 

4.2 Illustrative Cases 

The above survival calculations are presented as a guide to 
indicate which regimes are likely to be the most interesting, 
in terms of evaporation having an important role on 
planetary evaporation. In this section however we construct 
more accurate models that account for the evolution of the 
planet mass as well as the transition from X-ray driven 
to EUV driven evaporation. Therefore, we consider the 
two cases: that of a Jupiter mass planet and that of a 
Neptune mass planet, at orbital separations of 0.0125, 0.025, 
0.05 and 0.1 AU, and planetary densities of 0.1, 0.5 and 
1 g cm" 3 . Following the radial evolution of the planet is 
beyond the scope of this work, since this would require 
coupled modelling of the planetary interior, and thus we take 
the planets to evolve at constant density. 

Furthermore, should the planet transition from X-ray 
driven to EUV driven, we also check to see whether the 
EUV evaporation is in the 're combination limit' or th e 
'energy limit' using the results of iMurrav-Clav et al.l (|2009l ). 
Specifically, if the recombination-limited EUV flow gives a 
mass-loss rate above that of the energy-limited case, we 
switch to t he energy-limited formula for EUV evaporation 
provided bv lMurrav-Clav et all l|2009l 'l . In the energy-limited 
case the mass-loss rate will now scale linearly with the 
ionizing luminosity, rather than as the square root as in the 
recombination-limited case. 




10 100 

Planetary Mass [MJ 

Figure 10. Crude planetary survival lines at different separations 
compared to the observations. The points and lines are coded 
correspondingly: the solid blue line indicates the planetary 
survival line at 0.1 AU and blue diamonds the known planets 
with separations 0.05-0.1 AU, the survival line for 0.05 AU is red 
dashed and planets at 0.025-0.05 AU are red circles, the survival 
line at 0.025AU is black dotted and planets at 0.0125-0.025 AU 
are black squares, and the survival line for 0.0125AU is magenta 
dot-dashed while there are no planets at < 0.0125 AU. Planets 
above their respective survival line in the figure should have been 
destroyed according to the models. We also show lines of constant 
density as a guide at densities of 0.2, 1 and 5 g cm -3 . 

The results of these evolution calculations are shown in 
Fig. 1121 where the upper panels show planets with an initial 
mass of one Jupiter mass, and the lower panels show planets 
with an initial mass of one Neptune mass with separation 
decreasing from left to right. Note that planets with a density 
of 0.1 g cm -3 at a separation of 0.0125 AU are not shown 
as these planets would be undergoing dynamical Roche-lobe 
overflow rather than evaporation. In addition, a Jupiter mass 
planet with a density of 1 g cm" 3 at 0.1 AU does not 
experience hydrodynamic evaporation, and is not shown. 

These calculations show that Jupiter mass planets 
experience relatively little evaporation, losing no more than 
10 per cent of their initial mass, even at the closest 
separations. However, Neptune mass planets are prone to 
destruction at separations < 0.05 AU and can loose several 
tens of per cent of their initial mass at separations of 
0.05 — 0.1 AU. Thus we would not expect to see Neptune 
mass planets with low densities at small separations, a s 
indicated by the observations (e.g. [Hansen fc MurravjkoilT ). 
At late times we also notice a clear distinction in the driving 
mechanism as a function of separation. At separations < 
0.025 AU the evaporation of Jupiter mass planets appears 
to be X-ray driven while at larger separations it is EUV 
driven. This prediction could be compared to observations 
of the neutral content of the wind as a function of separation, 
though it must still be remembered that an X-ray driven flow 
can still have an EUV heated and fully ionized region above 
the sonic surface. 

4-2.1 Comparison with previous studies 

There have been several previous attempts to quantify the 
role of evaporation on the evolution of close in planets, 
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Figure 11. The X-ray luminosity at which a planetary wind will 
transition from being X-ray driven to EUV driven is shown for 
different separations. The points indicate Jupiter mass planets 
(unfilled) and Neptune mass planets (filled), with densities of 0.1 
(black squares), 0.3 (red diamonds) & 1 (blue stars) g cm" 3 . 



Lammer et al. (2009) u se the 'energy-limi ted' formalism for 
pure EUV evaporation, IPenz et al.l |200si ) use fixed heating 
effic iency models t hat in clude X-ray and EUV luminosity 
and I Jackson et aD (|2012T ) use the 'energy-limited' formalism 
to assess the role of pure X-ray evaporation. 

Lammer et al. (2009) estimated that pure EUV 
evaporation with an efficiency of 10% was not capable of 
completely destroying most close in planets, with Jupiter 
mass planets being more stable than lower mass planets, 
similar to what is found in this work. While the inclusion of 
X-ray mass-loss at earlier times does increase the total mass 
lost by close-in planets, it is not sufficient to make Jupiter 
mass planets unstable to complete evaporation. However, 
we do find lower mass (~ Neptune mass) planets are more 
unstable to evaporation than predicted by Lammer et al. 
(2009). 

IPenz et alj (|2008l ) considered the evolution of hot 
Jupiter's similar to HD209458b and found total mass-loss 
over Gyr's of evolution higher than found here, this is 
simila r to the results of certain cases in the I Jackson et al.l 
(|2012l ) models for pure X-ray evaporation. This a rises for two 
reaso ns: i) the fixed heating efficiencies ( 60% for lPenz et al.l 
120081 and 25% for Ijackson et all 120121 ) are rather large 
when compared to the actual values obtained in certain 
regimes (see Section 15,2.20 : ii) the X-ray and EUV field 
were considered to both be contributing to the heating at 
all times, whereas we have shown there is a switch from 
X-ray driven to EUV driven at late times which results in 
reduced mass-loss rates. Thus we conclude that estimates 
of total mass-loss at levels around 30% for hot Jupiter's at 
separations of 0.01-0.02AU are unrealistically high. 

The use of X-ray evaporation at earlier times which 
switches to EUV evaporation at late times thus represents 
an intermediate between pure EUV evolution models, and 



models which include a contribution from the X-ray at all 
times. 



4-2.2 Limitations 

The major limitation of these evolution calculations is that 
we do not follow the evolution of the interior structure of 
the planet and hence its radius. We make the assumption 
the planet evolves at constant density, motivated by the 
fact that the spread in the density of observed planets is 
comparatively small relative to the orders of magnitude 
changes in mass. It is, however, worth noting that a planet 
whose density decreases with time (as is the case for many 
mode ls of planets undergoing mass loss, e.g. iBaraffe et al.l 
|2004| ) will be: i) prone to higher mass-loss rates for a given 
mass; and ii) experience a transition from X-ray driven to 
EUV driven at later times. These effects will result in more 
rapid planetary evaporation than our calculations here. On 
the other hand if the density of the planet increases with 
time, as may happen if a large fraction of the planet mass is 
comprised of a rocky core, the converse will be true and the 
planet may be able stabilize itself ag ainst evaporation. 

In addition. lBaraffe et al.l (|2004 ) point out that a critical 
point may occur when a planet's evaporation timescale 
becomes shorter than the thermal timescale. In this case 
the planet is no longer able to shrink its radius to adjust 
to the mass-loss, and the planets evaporation is accelerated 
(lower mass with same radius results in higher mass-loss) 
and can lead to runaway mass-loss. This effect, particularly 
combined with a general tendency for density to decrease 
with decreasing mass, may make planets more prone to 
destruction compared to our simple calculations. 



5 DISCUSSION 

In the previous sections we have performed basic 
hydrodynamic calculations to investigate the evaporation 
of close-in planets, and the various limits it may exist in. 
Qualitatively we have shown that most close-in planets (a < 
0.1 AU) will experience a period of hydrodynamic mass-loss 
driven by X-rays at early times, followed by EUV driven at 
late times, as the high-energy luminosity drops off with time. 
Furthermore, due to the steep nature of the X-ray heated 
sub-sonic flow, the atmospheric structure is unimportant in 
determining the flow structure, although as a consequence 
the flow structure will be sensitive to the planetary radius. 
By considering the combination of EUV and X-ray heating, 
we have shown that X-ray evaporation dominates when the 
sonic surface occurs in the X-ray heated flow, typically for 
high mass-planets with large radii and a high X-ray flux. 
When the EUV evaporation dominates on the other hand, 
the X-ray portion of the flow is entirely sub-sonic, and EUV 
domination typically occurs for low mass-planets with small 
radii and low X-ray flux. 



5.1 Hydrodynamic or Ballistic Escape? 

One of the questions we set out to answer in this work 
is whether, in general, the evaporation of close-in planets 
occurs through a hydrodynamic process, or through a 
ballistic process (e.g. Jean's escape). It is clear from this 
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Figure 12. The mass evolution of close-in planets for initially Jupiter mass planets (upper panels) and Neptune mass planets (lower 
planets) at densities of 0.1, 0.3 & 1 g cm -3 . Planets of lower density lose mass faster. Note that planets of density 0.1 g cm -3 at 0.0125 AU 
and Jupiter mass planets of density 1 g cm -3 at 0.1 AU are not shown as the flow is not in the hydrodynamic limit. From left to right 
the panels indicate decreasing separations from 0.1, 0.05, 0.025 to 0.0125 AU. Dashed black lines indicate the time when the evaporation 
is dominated by X-rays and solid red lines when the evaporation is dominated by the EUV. 



work that the dominant mass-loss process for close-in planets 
(a < 0.1 AU) is through a hydrodynamic wind. However 
above a certain mass dense planets are too massive to 
loose gas hydrodynamically, either through X-ray or EUV 
evaporation. Fig.[6]and[9]indicate that planets more massive 
than 1 Jupiter mass at densities greater than 1 g cm -3 
have atmospheres which are too tightly bound to drive a 
hydrodynamic wind at separations ~1 AU. At very close 
orbital separations however, < 0.05 AU, even planets of 
a few Jupiter masses are able to undergo hydrodynamic 
evaporation. As the X-ray luminosity drops, so does the 
region which is prone to hydrodynamic evaporation, as 
clearly indicated by Figure[6] Given the much higher collision 
cross-section of gas particles in fully ionized (and therefore 
EUV heated) gas, planets switch to EUV driven evaporation 
and remain EUV driven before they become unable to 
evaporate hydrodynamically for all Jupiter and Neptune 
mass planets considered in Section [4] 



5.2 Limitations of the models 

These models have given valuable insight into which 
processes dominate the evaporation in different regions of 
parameter space. However, it is important to highlight 
the limitations of the models and in particular discuss 
the two major assumptions of this work: adopting a solar 
composition for the X-ray heating (although we do account 
for an overall metallicity scaling) and assuming radiative 
equilibrium. 



5.2.1 Role of atmospheric composition 

Perhaps one of the major uncertainties when discussing 
X-ray flows is the role of composition, since X-ray heating 
and cooling is done by metals. Although calculations of 
atmos p heric composition ar e begi n ning to be per f ormed 
dYelld 12004 iGarcfa Munozl 120071 ; iKoskinen etafl 12001 
120101 ). as a starting point we have assumed the metals have 
abundance ratios comparable to solar values, and included 
a basic metallicity s caling derived from pre vious work 
on evaporating flows (jErcolano fc Clarke! I201CI ) . Thus our 
models are unlikely to be appropriate for considering very 
exotic, metal rich atmospheres or where the compositions 
are significantly different from solar. The discrepancy may 
not be as bad as one would originally expect though, since 
as discussed in Section [2] X-ray photo-electron heating, and 
line cooling, are typically done by the same metal species, 
with O xygen and Carbon t ypically being the dominant 
species (|Ercolano et al.ll200ct ). As Oxygen and Carbon has 
been detected in evaporating atmospheres of close-in planets 
(e.g. IVidal-Madiar et ail |2004 | Ben-Jaffel fc Sona Hosseinil 
|201C( ). this is reassures us that the models presented in this 
work represent a good starting point. Nonetheless future 
work is certainly needed to determine the role atmospheric 
composition has on evaporation, along the lines of coupling 
the photo-chemistry of important species e.g. Hydrogen, 
Helium, Carbon and Oxygen to the hydrodynamic models. 



5.2.2 Energy efficiency 

Estimates of hydrodynamic evaporation have previously 
been calculated through so called 'energy-limited' 
evaporation where most of the received energy, assumed to 
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Figure 13. The 'efficiency' of the X-ray wind calculated by 
comparing the energy input received by the wind out to 
the sonic surface, to the mechanical luminosity of the wind 
[4L moc i 1 a 2 /(Lx^s)l ■ Here the calculation is performed for the 
X-ray flow of Fig. \E\ though note that as -t/ mec h °c Lx/a? this 
efficiency is invariant under changes to Lx and a provided the 
change in R B is small. 

be deposited at the planet surface, can be taken to be lost 
in hydrodynamic expansion as PdV work (e.g. Watson et 
al. 1981). Howev er, there is no a pr i ori re ason to expect 
this, and in-fact iMurrav-Clav et"al] (|2009l ) determined a 
transition between energy-limited EUV evaporation at low 
EUV fluxes, and the case where the received energy is 
predominately re-emitted as cooling radiation at high EUV 
fluxes. Furthermore, comparisons of the heating and cooling 
timescales with the dynamical timescales of X-ray heated 
gas in this ionization parameter range indicate that the 
X-ray gas is close to radiative equilibrium (Owen et al. 
2010). Thus similarly for X-ray winds, the energy emitted 
as cooling radiation is considerably larger than that which 
goes into PdV work, suggesting much lower efficiencies than 
typically assumed using the energy-limited method. This 
may not ultimately result in much lower evaporation rates 
however, as the X-ray wind is free to absorb energy out to 
its sonic surface. 

In this work we have assumed that the X-ray gas 
is in radiative equilibrium and neglected PdV cooling 
based on comparisons of the dynamical and thermal 
timescales. However, it is prudent to check this assumption 
by comparing the mechanical luminosity (L mec h = rn w X 
Ae, where e is the total specific energy of the gas) of 
our X-ray evaporative winds to the energy input rate 
and calculate the 'efficiency' of our flows. An efficiency 
value > 1 would indicate that our assumption of radiative 
equilibrium is poor and that cooling due to PdV" work 
is important. In Fig. [13] we plot the wind 'efficiency' 
(4L mec ha 2 /-L;sc Rs) for the flow shown in Fig. [S] In general 
this efficiency is low (< 0.1) for planets with Neptune 
masses and above, indicating our assumption of radiative 



equilibrium is a good one, though we note that the efficiency 
increases to slightly higher values, ~ 0.15, for high density, 
earth-like, planets. Jupiter mass planets drive the least 
efficient winds, with efficiencies of ~ 0.01. This increase 
in efficiency with decreasing mass is most likely caused 
by the decrease in the dynamical timescale for planets 
with lower mass, whereas the thermal timescale remains 
roughly constant. Additionally given that the mass- loss 
rate (and hence mechanical luminosity) scales linearly 
with X-ray flux, this indicates that X-ray evaporation is 
close to radiative-equilibrium at all X-ray luminosities and 
separations and that unlike the EUV case there is no 
transition from radiation-limited to energy-limited. Thus, 
our calculations assuming radiative-equilibrium are accurate 
for the parameter space considered, and X-ray evaporation 
should not be thought of as 'energy-limited'. Furthermore, 
while we point out that for sensible choices of the efficiency 
value, the 'energy-limited' formalism can be used to recover 
correct values of the mass-loss rate, and give similar scalings, 
we caution that this efficiency value is certainly not fixed over 
orders of magnitude in both planetary mass and radius, and 
it cannot take into account how the increased area due to 
the position of the sonic surface varies with the parameters. 
Therefore, we find that Eq. fJJ is of limited usefulness 
and would advocate using more accurate estimates of the 
mass-loss rates. Furthe rmore, comparing thes e effici encies 
with the calculations of ICecchi-Pestellini et all HQOS), who 
calculate various heating efficiencies for fixed ionization 
fractions (10 _1 — 10 -6 ) and find efficiencies between 80% 
and 2% respectively, we see that in all cases the bulk of 
the flow is at a fairly low ionization level. Calculating the 
ionization fraction directly from the ionization parameter 
in our calculations we find ionization fractions > 10~ 2 in 
the transonic region of the flow, in agreement with the 
static calculations of ICecchi-Pestellini et all |2006l ). further 
validating the ionization parameter approach used in this 
work. 

5.2.3 Multi-dimensional problem 

Planetary evaporati on is ultimately a 3D problem (e.g. 
IStone fc Progall2009l ). and comparisons with observed line 
profiles will require knowl edge of the 3D structure of the flow 
(|Murrav-Clav et al.ll2009 l). and the role of day vs night-side 
heating. The models in this paper represent formal upper 
limits for the mass-loss rates, as we have assumed spherical 
outflow over the full 4ir angular area of the planet. Future 3D 
radiation-hydrodynamical models of planetary evaporation 
will be able to shed light on the 3D structure of these flows 
(Owen et al. in prep), which will allow comparisons with 
observed line profiles, as lines sensitive to the flow outside 
the Roche radiu s will certainly be aff ected the 3D nature 
of the problem |Koskinen et al.l |2010| ). Furthermore, they 
will be able to locally test some of the approximations (e.g. 
radiative equilibrium), as well as investigate the interaction 
between EUV and X-ray heating in the 3D setting. In this 
work we have assumed that the transition between X-ray 
driven and EUV driven occurs at a fixed point in time 
across the entire planetary wind, however, it is possible 
that parts of the planet could be undergoing EUV driven 
evaporation while other parts are X-ray driven. Furthermore, 
a full 3D model will allow the assessment of whether 
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non-axisymetric mass-loss can play a role in the orbital 
evolution of the planet, particularly for those undergoing 
significant frac tional mass- l oss ( Neptune mass planets) as 
investigated bv lBoue et ail (|2012l ). 



6 CONCLUSIONS 

In this work we have investigated the role of both X-ray and 
EUV heating in evaporating a close-in planet's atmosphere, 
using ID hydrodynamic models. We have considered how the 
mass-loss rates, and the driving mechanism, vary within the 
observable parameter space including: planet mass, planet 
density, separation, stellar mass, metallicity of the planet 
atmosphere and X-ray/EUV luminosity. Then under the 
assumption that the EUV and X-ray luminosities follow an 
identical evolution in time, with the same initial luminosities, 
we have calculated the mass evolution of evaporating planets 
that evolve at constant density. Our main conclusions are as 
follows: 

(i) We identify two separate cases of planetary 
evaporation: X-ray driven, where the X-ray flow undergoes 
a sonic transition, and EUV driven, where the X-ray 
heated portion of the flow remains sub-sonic. In general 
X-ray evaporation occurs at high X-ray luminosities, low 
planetary densities, high planetary masses and small 
separations. Whereas EUV evaporation dominates at low 
X-ray luminosities, high planetary densities, low planetary 
masses and large separations. 

(ii) We find that at separations of <0.1 AU, most planets 
will be evaporating hydrodynamically rather than through a 
ballistic mass-loss process. However, at separations greater 
than ~ 0.1 — 0.5 AU dense, Jupiter mass planets, may not 
be able to evaporate hydrodynamically. 

(iii) In the case of X-ray driven evaporation, the flow is 
close to radiative equilibrium and energy loss from the flow is 
dominated by line cooling. The mass-loss rates scale linearly 
with X-ray luminosity and inversely with the separation 
squared (m oc Lx/a 2 ), but no exact scalings with planetary 
mass and planetary radii exist. 

(iv) We find that as the high-energy luminosity falls over 
time, the evaporative flow may undergo a transition from 
X-ray driven at early times to EUV driven at late times. 
This transition occurs at lower X-ray luminosities (and hence 
later times) for planets with at smaller separations. 

(v) Considering the evolution of Jupiter and Neptune 
mass planets, we find it is unlikely that an initially Jupiter 
mass planet can be completely evaporated within Gyr 
timescales, only losing a few percent of its original mass. 
However, we find that Neptune mass planets are much more 
susceptible to complete evaporation. Within 3 Gyr Neptune 
mass planets at densities of <0.3 g cm~ 3 are completely 
evaporated at a separation of 0.025 AU, while Neptune mass 
planets at densities of <1 g cm -3 are completely evaporated 
at a separation of 0.0125 AU. 

(vi) Using our evolutionary models we infer that at late 
times (~Gyrs), planets at separations <0.025AU will still 
be undergoing X-ray driven evaporation, whereas at larger 
separations the evaporation will generally be in the EUV 
driven regime. 



ACKNOWLEDGMENTS 

We thank the anonymous reviewer for comments which 
helped improve the paper. AJ is supported by an STFC 
Postgraduate Studentship. We thank Barbara Ercolano, 
Cathie Clarke, Norm Murray and Yanqin Wu for productive 
discussions. JO is grateful to hospitality from the IoA, 
Cambridge where this work began. 



APPENDIX A: GENERAL X-RAY FLOW 
SOLUTION 

In the case of a general X-ray heated flow, where the 
gas temperature is specified by a monotonic function of 
the ionization parameter (T = /(C)), we can follow a 
similar analysis to that given in Section 3, which gives 
an equation for the critical sub-to-super-sonic transition 
provided dlog//dlog£ < 1 (a < 1 as in the illustrative 
case shown in Section 3.1) as 



dlog/ 
dlogC 



f=f< 



2 dr 
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where Eq. IA2I can be solved for the ionization parameter 
at the sonic surface (£„). Using these conditions to specify 
the ionization parameter at a chosen sonic surface, as 
before we can obtain a transcendental function for the 
solution of £(r), by applying mass conservation and using 
the Bernoulli potential. Thus the stream-function for the 
streamline connecting the planet and star is 
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where H(R S ) is the Bernoulli constant for a flow with a sonic 
transition at that given radius or 
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Eq, IA3l can now be solved for a family of solutions dependant 
on R s for a specified T = /(£). 
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